Goto

Collaborating Authors

 use mcmcglmm


Using R: Animal model with simulated data

#artificialintelligence

Last week's post just happened to use MCMCglmm as an example of an R package that can get confused by tibble-style data frames. To make that example, I simulated some pedigree and trait data. Just for fun, let's look at the simulation code, and use MCMCglmm and AnimalINLA to get heritability estimates. Now that we have some data, we can fit a quantitative genetic pedigree model ("animal model") to estimate genetic parameters. We're going to try two methods to fit it: Markov Chain Monte Carlo and (the unfortunately named) Integrated Nested Laplace Approximation.